Violation of self-similarity in the expansion of a ID Bose gas 
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The expansion of a ID Bose gas is investigated employing the Lieb-Liniger equation of state 
within the local density approximation. We show that during the expansion the density profile of 
the gas does not follow a self-similar solution, as one would expect from a simple scaling Ansatz. 
We carry out a variational calculation, which recovers the numerical results for the expansion, the 
equilibrium properties of the density profile, and the frequency of the lowest compressional mode. 
The variational approach allows for the analysis of the expansion in all interaction regimes between 
the mean field and the Tonks-Girardeau limits, and in particular shows the range of parameters for 
which the expansion violates self-similarity. 



I. INTRODUCTION 

The experimental achievement of Bosc-Einstcin con- 
densation (BEC) Q has aroused a large interest in the 
physics of ultracold atomic gases. Among the topics re- 
lated to this field, the physics of low dimensional atomic 
gases has recently attracted significant attention. The 
development of the trapping techniques has allowed for 
the realization of very anisotropic geometries, where the 
confinement is so strong in one or two dimensions, that 
at low temperatures the transversal motion is "frozen" , 
and does not contribute to the dynamics of the system. 
In this way two- 3, 4, 5] and onc-dimcnsional 
systems have been accomplished. Low-dimensional gases 
present significantly different properties compared to the 
three-dimensional ones. A remarkable example is pro- 
vided by the existence of quasi-condensation 8..9. .1CL.lli] . 
whose effects have been recently observed experimentally 

m 

During the last years, the ID Bose gases have been 
the subject of growing interest, in particular the limit 
of impenetrable bosons fT3!| , which behave to a large ex- 
tent as a noninteracting Fermi system, acquiring some 
remarkable properties. The conditions for the experi- 
mental realization of strongly correlated ID gases are 
rather restrictive |^ Q], since a large radial compres- 
sion, a sufficiently small density, and eventually a large 
scattering length are needed. Fortunately, recent exper- 
imental developments have opened perspectives in this 
sense. Especially interesting is the possibility to modify 
at will the interatomic interactions by means of Feshbach 
resonances , and the cap ability of loading an atomic 
gas in an optical lattice 

From the theoretical side, the physics of ID Bose gases 
was first investigated by Girardeau 0, who consid- 
ered the limit of impenetrable bosons, also called Tonks- 
Girardeau (TG) gas, pointing out a non trivial relation 
with the physics of ideal Fermi gases. This analysis was 
later extended by Lieb and Liniger "2l], who solved ana- 
lytically the problem for any regime of interactions, using 
Bethe Ansatz. Yang and Yang ||22| extended the anal- 
ysis including finite temperature effects. Recently, the 



experimental accessibility of trapped gases, have encour- 
aged the investigation of the harmonically trapped case. 
The Bose-Fermi (BF) mapping has been empl oyed to the 
case of an inhomogeneous gas in the TG limit '23^. How- 
ever, there is unfortunately, to the best of our knowl- 
edge, no exact solution for arbitrary interaction strength 
in the case of trapped gases. The problem of the equi- 
librium of a trapped gas can be analyzed using a local 
density approximation and employing the Lieb-Liniger 
(LL) equation of state locally to evaluate the equilibrium 
density profiles ^17]. A similar formalism has been re- 
cently employed to analyze the collective oscillations in 
the pres ence of harmonic trapping 0. Both Refs. 
and [ig have shown the occurrence of a continuous tran- 
sition from the mean field (MF) regime to the TG one as 
the intensity of the interaction is varied. Recently, Gan- 
gardt and Shlyapnikov have discussed the stability 
and phase coherence of ID trapped Bose gases. These au- 
thors have analyzed the local correlation properties and 
found that inelastic decay processes, such as three body 
recombination, are suppressed in the TG regime, and in- 
termediate regimes between MF and TG. This fact opens 
promising perspectives towards the accomplishment of 
strongly interacting ID Bose gases with large number of 
particles. This analysis have been very recently extended 
to the case of finite temperatures 0| . 

The expansion of a one-dimensional Bose gas in a guide 
was analyzed in Ref. pol |. by means of a hydrodynamic 
approach based on the local LL model. The expansion 
dynamics was shown to be different for different inter- 
action strengths, and its analysis could be employed to 
discern between the TG and MF regimes. In particular, 
the self-similar solution is violated. 

In this paper, we extend the analysis of Ref. '2^1 by 
introducing a variational approach, which permits us to 
study the asymptotic regime at large expansion times. 
This method is shown to be in excellent agreement with 
previous direct numerical simulations, and additio nally 
permits us to recover the results of Refs. 17] and [l8l |. 
More importantly, our variational approach allows us to 
determine the regime of parameters for which the self- 
similarity of the expansion is violated. 

The paper is organized as follows. In Sec. ^1 we in- 
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troduce the local LL model which we employ to analyze 
the expansion dynamics. In Sec. 11111 we briefly discuss 
the numerical results obtained in Ref. j23|. In Sec. II VI 
we present a variational approach which allows us to in- 
vestigate in detail the expansion dynamics for arbitrary 
regimes of parameters. Finally we conclude in Sec. 



of equations [2l| 



Hi) = 7 / 9{xh)dx. 



(6) 



II. LOCAL LIEB-LINIGER MODEL 

We analyze in the following a dilute gas of N bosons 
confined in a very elongated harmonic trap with radial 
and axial frequencies ujp and lo^ (up ^ lOz)- We assume 
that the transversal confinement is strong enough so that 
the interaction energy per particle is smaller than the 
zero-point energy hujp of the transversal trap. In this 
way, the transversal dynamics is effectively "frozen" and 
the system can be considered as dynamically ID. In this 
section we briefly review the formalism introduced in Ref. 

m- 

We assume that the interparticle interaction can be ap- 
proximated by a delta function pseudopotential. There- 
fore the Hamiltonian that describes the physics of the ID 
gas becomes 



N 



H 



ID 
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(1) 



where 



^ID 



2m ^ 



N-l N 

3iD ^ ^ 5{z,-Zj) (2) 



is the homogeneous Hamiltonian in absence of the har- 
monic trap, m is the atomic mass and (/id = —2h^/maijy. 
The scattering problem under one-dimensional con- 
straints was analyzed in detail by Olshanii [T^ . and 
it is characterized by the one-dimensional scattering 
length oiD — [—a?p/2a)[l — C{a/ap)\, with a the three- 
dimensional scattering length, Cp = ^^hjraiijp the oscil- 
lator length in the radial direction, and C = 1.4603 .... 
As shown by Lieb and Liniger |2ll |. the homogeneous 
Hamiltonian g? T-, ca n be diagonalized exactly by means 
of Bethe Ansatz ^5]. In the thermodynamic limit, a ID 
gas at zero temperature with a given linear density n, is 
characterized by the energy per particle 



where 7 — 2/n|aiD|. The function 5(7) fulfills 



e(7) 



A3(7) 



g{xYi) x'^dx, 



(3) 



(4) 



We assume next that at each point z the gas is in local 
equilibrium, with local energy per particle provided by 
Eq. (O. Then, one can obtain the corresponding hydro- 
dynamic equations for the density and the atomic veloc- 
ity 



9 d , , 
ot oz 



(7) 



where 



liie{n) = ( 1 + n— ] e(n) 



(9) 



is the Gibbs free energy per particle. 

Note that for the case of n|ai£i| — > cx), one obtains 
Mie {n) = g\Dn;j:BtT:'ieYm.g the ID Gross-Pitaevskii equa- 
tion (GPE) 12(1, whereas for the case Tijaid ~* 0) ^^^^ 
gets niein) = 7r^/i^n^/2m, and the equation of Ref. 
|27| is recovered. The system has only one control pa- 
rameter 

B in m, namely A = A^jainp / a1 , where 
Oz = ^/hJmjjTz is the harmonic oscillator length in the 
z direction. The regime A 3> 1 corresponds to the MF 
limit, in which the stationary-state density profile has a 
parabolic form. On the other hand, the regime A <C 1 
corresponds to the TG regime, which is characterized 
by a stationary-state density profile with the form of a 
square root of a parabola. 



III. NUMERICAL RESULTS 

In Ref. Il^l, equations © , I© , © , I© were employed 
to simulate numerically the expansion of a ID gas in the 
framework of the hydrodynamic formalism. The expan- 
sion follows the sudden removal of the axial confinement, 
while the radial one is kept fixed. In particular, it was 
observed that during the expansion the density profile is 
well described by the expression 



n(z,i) = 1 - 



(10) 



where g(x\^^ and A (7) are the solutions of the LL system 



where nm{t) provides the appropriate normalization, 
R{t) is the radius of the cloud, and the exponent s(t) 
takes the value s(0) = 1 for an initial MF gas. The func- 
tion s{t) decreases monotonically in time, approaching 
an asymptotic value (see Fig. QJ. Therefore, contrary 
to the expansion dynamics for a BEG [2^ 123 . l30| , the 
self-similarity of the density profile is violated. 
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At this point we discuss the physics behind this vio- 
lation of the self-similarity. If the local chemical poten- 
tial presents a fixed power law dependence on the den- 
sity, (X n^, it is easy to show from the hydrodynamic 
equations ^ and ((SJ, that there exists a self-similar so- 
lution of the form n — (no/6)(l — {z /hRYY^^ ^ where 
h = Lol/h^^^ . For the particular case of the TG gas, the 
scaling law can be also obtained from the exact BF map- 
ping [2^. However, since /i/e is obtained from the LL 
equations, the dependence of /i;e on n is quadratic for a 
low density and linear for a large one. Therefore, /i;e does 
not fulfill a fixed power law dependence during the ex- 
pansion, and the self-similarity is violated. In particular, 
as the expansion proceeds the whole system approaches 
the low density regime, and consequently the exponent s 
decreases monotonically. In the next section, we analyze 
in more detail this effect. 



assume the following Ansatz for the density 




timeCl/co^) 



FIG. 1: Time evolution of the exponent s{t) for A = 0.43, 
ujp = 2-n X 20kHz and N = 200 atoms {uo^ = 2ti x 1.8Hz at 
t = Qi). Our variational result (dashed line) shows a very good 
agreement with the results obtained from the direct resolution 
of Eqs. and 



IV. VARIATIONAL CALCULATION 



C{s) 



1 - 



(12) 



where h and s are time dependent variables, R is the 
initial Thomas-Fermi radius, and C(s) is related to the 
normalization to the total number of particles. For the 
<j) field we consider the following form: 



1 



1 



-az 



-Pz" 



(13) 



where a and (3 are time dependent parameters. We stress 
at this point, that in the analysis of the self-similar ex- 
pansion of a BEG |33, a quadratic ansatz (in 
z) for the field provides the exact solution. However, 
for the problem under consideration in this paper, it is 
necessary to include higher order terms to account for 
the violation of the self-similarity. We have checked that 
terms of higher order than introduces only small cor- 
rections, and therefore we reduce to the form of Eq. I|13|l . 




ln(A) 



FIG. 2: Equilibrium values (at t = Q before opening the trap) 
of the exponent s, as a function of the parameter A. The 
dashed lines denote the MF limit, so = 1, and the TG one, 
so = 1/2. 



In this section, we complete our understanding of the 
expansion of a one-dimensional Bose gas in a guide by 
means of a variational Ansatz using a Lagrangian for- 
malism. The Lagrangian density for the system is of the 
form 



90 1 

L = —mn— mn 

dt 2 



— ^ ) miJlz^n — ein^n (11) 

dz 2 " \ ) \ ) 



where the velocity field is defined as w = d(j)/dz. The 
equations of motion are obtained from the functional 
derivation of the action A = j Cdtdz: SA/d(j) — (conti- 
nuity equation), SA/6n = (which after partial deriva- 
tion with respect to z provides the Euler equation) . From 
the numerical results we have observed that the density 
is at any time well described by Eq. H1U|) . Therefore, we 
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FIG. 3: Frequency of the lowest compressional mode as a 
function of the parameter A. 
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We are interested in the dynamics of the parameters 
b and s, related to the size and the shape of the cloud, 
respectively. Integrating the Lagrangian density in 2, 
L = J Cdz, one finds a Lagrangian for the above men- 
tioned parameters: 



mNR^ 




(3b^ 



2s + 3 2 (2s + 5)(2s + 3) 
a(3b* 



2s + 3 



(2s + 5)(2s 



3) 



(2s + 7) (2s 
b^Lo: 



2s + 3 



5)(2s + 3) 
dzne{n). 



We perform a gauge transformation jsl 

d 



L{t, q, q) L{t, q, q) + ^.9(^, q), 



(14) 



(15) 



where 



(ib'^ 



2 [ 2s + 3 2 (2s + 5)(2s + 3) 

■ (16) 

The resulting Lagrangian is of the form L — 
L(a, P,b,b^ s, s). Imposing the conservation laws 
dL/da = dL/d(3 = 0, we obtain a Lagrangian depending 
only on the two relevant parameters s and 5, of the form 
L = K — V , where 



V 

EiD 

K 

EiD 



B{b,s) 
/o(s) 
(Mnb^ 



2l2 



A'b 



1 



[77o/o(so)]2 2s + 3 
+ 2Mi2bs + M22S' 



[?7o/o(so)]2(2s + 3) 



(17) 



(18) 



where i^iD = fi,^/2m|ai£)p is the typical energy asso- 
ciated with the interatomic interactions. In Eq. (|17|l . 
we use the function B{b, s) = J dy{l — y'^Y e{n{y)) / Eijj , 
where we integrate over the rescaled axial coordinate 
y — z/Rb. In Eqs. (|17|l and (|18|l . we define the dimen- 
sionless central density ryo = '^oI'^I-dL where no is the 
initial central density, and the parameter sq — s{t — Q). 
We have additionally employed the auxiliary functions 
fn{s) — J y"(l — y'^ydy, and the coefficients 



Mil 
M12 

Af22 



= 1, 



_b_ 

27+3 
62(121 



186s + 96s2 + 16s3) 



4(s-f l)(2s-f3)2(2s -1-5)2 



(19) 
(20) 

(21) 



From the Lagrangian L, one obtains the corresponding 
Euler-Lagrange equations for the parameters b and s. In 
order to find the initial conditions 770 and sq for the ex- 
pansion, we have numerically minimized the potential V 



in the presence of the harmonic trap for different values 
of A, assuming b —\ (see Fig.|21). When \, so tends 
to 1, as expected for the MF case. On the contrary, when 
A ^ 1, So tends to 1/2 (TG proffle). As expected from 
Ref. [l3, for A < 1, 770 oc A^/'^, whereas for A > 1, the 
MF dependence 770 oc A^/'^ is recovered. 

Our variational approach allows us to calculate 
the lowest compressional mode, offering an alternative 
method as the one discussed in Ref. ^8| . Expanding the 
potential V around the equilibrium solution up to second 
order in b (see Fig. O, and neglecting for small oscilla- 
tions the time-dependence of s, we obtain 



1 HfMY d'B 



2^2 /2(so) 



952 



(22) 



6=1 



Our results show a continuous transition from the MF 
value. Urn = VSuJz, to the TG one, LUm = '2llIz, in excellent 
agreement with the results obtained by means of a sum 
rule formalism [Tsf . 
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FIG. 4: Expansion velocity. Asymptotic value of & as a func- 
tion of the parameter A. 



ln(A) 



FIG. 5: Violation of the self-similarity. Value of ^ = Soo/so 
as a function of the parameter A. The value ^ = 1 denotes 
self-similarity. 

From the corresponding Euler-Lagrange equations, we 
have obtained the dynamics of b{t) and s{t). We have 
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checked in all our calculations that the energy and num- 
ber of particles remain a constant of motion. We have 
compared the variational results with our simulations 
based on the exact resolution of the hydrodynamic equa- 
tions j22|], obtaining an excellent agreement (see Fig.^. 

We have analyzed the asymptotic value, 600, for differ- 
ent values of A. It is easy to obtain, that for a power law 
dependence of the local chemical potential fii^ oc , the 
derivative of the scaling parameter asymptotically ap- 
proaches a value — \/^T^i^z- Therefore a continuous 
transition from 600 — (TG) to hoa = ^J^^z (MF) is ex- 
pected. We recover this dependence from our variational 
calculations (see Fig. 0J). 

We have analyzed the behavior of s during the expan- 
sion dynamics for different values of A. In particular, 
we have defined the asymptotic ratio ^ = Soo/sq, with 
Soo — s{t 00) (see Fig. El. Deeply in the TG regime 
(A < 1) or in the MF one {A > 1), ^ ~ 1, i.e. in those ex- 
treme regimes, the expansion is well-described by a self- 
similar solution. However, for intermediate values, ^ < 1, 
i.e. the expansion is not self-similar. The self-similarity 
is maximally violated in the vicinity of ^ = 1, although ^ 
departs significantly from 1 for a range 0.01 ^ A < 100. 
The behavior at large A can be understood as follows. If 
the gas is at i = deeply in the MF regime, the change of 
the functional dependence of nie with the density occurs 
at very long expansion times, when the initial interaction 
energy of the gas has been fully transferred into kinetic 
energy. Therefore, for large values of A the self-similarity 
is recovered. 



V. CONCLUSIONS 

In this paper, we have extended the analysis of Ref. [2^ 
on the expansion dynamics of a one-dimensional Bose gas 



in a guide. We have shown that the expansion violates 
under certain conditions the self-similarity, and in this 
sense differs significantly from the expansion dynamics of 
a BEG. We have shown that the problem can be solved 
by employing the hydrodynamic approach, and the local 
Lieb-Liniger model. We have developed a variational ap- 
proach based on a Lagrangian formalism to study the ex- 
pansion for any regime of parameters. We have identified 
the possible physical situations at which self-similarity is 
violated. This should occur in a rather wide range of 
parameters (0.01 A < 100). The particular proper- 
ties of the expansion of a gas in the strongly-interacting 
regime could therefore be employed to discern between 
mean-field and strongly-interacting regimes. In addition, 
the asymptotic behavior of the expanded cloud could be 
employed to discriminate between different initial inter- 
action regimes of the system. 

Our discussion has been restricted to the analysis of the 
density properties. In fact the present formalism cannot 
describe the dynamics of the coherence in the system, i.e. 
we are limited to the diagonal terms of the corresponding 
single-particle density matrix. The description of the non 
diagonal terms lies beyond the scope of this paper, and 
requires other techniques of analysis [Tglls^l . 
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